Transition to Instability in a Kicked Bose-Einstein Condensate 
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A periodically kicked ring of a Bose-Einstein condensate is considered as a nonlinear generalization 
of the quantum kicked rotor. For weak interactions between atoms, periodic motion (anti-resonance) 
becomes quasiperiodic (quantum beating) but remains stable. There exists a critical strength of 
interactions beyond which quasiperiodic motion becomes chaotic, resulting in an instability of the 
condensate manifested by exponential growth in the number of noncondensed atoms. Similar behav- 
ior is observed for dynamically localized states (essentially quasiperiodic motions), where stability 
remains for weak interactions but is destroyed by strong interactions. 
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The classical kicked rotor is a textbook paradigm for 
dynamical chaos . The quantum kicked rotor has play 
an equally important role for the study of quantum chaos, 
for which a wide range of effects have been predicted |2( 
and observed in experiments Q. In recent years, the 
realization of Bose-Einstein condensation (BEC) of di- 
lute gases has opened new opportunities for studying 
dynamical systems in the presence of many-body inter- 
actions. A natural question to ask is how the physics of 
the quantum kicked rotor is modified by interactions 

In the classical regime, chaotic motion leads to diffu- 
sive growth in the kinetic energy. In quantum dynamics, 
where chaos is no longer possible because of the linearity 
of the Schrodinger equation, the motion becomes periodic 
(anti-resonance), quasiperiodic (dynamical localization), 
or resonant (quantum resonance) -"- n the mean field 

approximation, many-body interactions in BEC are rep- 
resented by adding a nonlinear term in the Schrodinger 
equation |8| • This nonlinearity makes it possible to bring 
chaos back into the system, leading to instability (in the 
sense of exponential sensitivity to initial conditions) of 
the condensate wave function |9( . The onset of instability 
of the condensate can cause rapid proliferation of ther- 
mal particles |l0| that can be observed in experiments 
[llj| . It is therefore important to understand the route 
to chaos with increasing interactions. This problem has 
recently been studied for the kicked BEC in a harmonic 
oscillator [l2j |. 

In this Letter, we investigate the quantum dynamics 
of a BEC with repulsive interaction that is confined on 
a ring and kicked periodically. This system is a nonlin- 
ear generalization of the quantum kicked rotor. From 
the point of view of dynamical theory, the kicked rotor 
is more generic than the kicked harmonic oscillator, be- 
cause it is a typical low dimensional system that obeys 
the KAM theorem [j^. It is very interesting to un- 
derstand how both quantum mechanics and mean field 
interaction affect the dynamics of such a generic sys- 
tem. We will focus most of our attention on the case 
of anti-resonance because it is the simplest periodic mo- 



tion. Here we find, with both analytic and numerical 
calculations, that weak interactions will make the motion 
quasiperiodic in the form of quantum beating. For strong 
interactions, quasiperiodic motions are destroyed and we 
observe a transition to instability of the BEC that is also 
characterized by an exponential growth in the number 
of noncondensed atoms. Universal critical behavior for 
the transition is found. We have also studied nonlin- 
ear effect on dynamically localized states that may be 
regarded as quasiperiodic. Similar results are obtained 
in that localization remians for sufficiently weak interac- 
tions but become unstable beyond a critical strength of 
interactions. 

Consider condensed atoms confined in a toroidal trap 
of radius R and thickness r, where r <C R so that lat- 
eral motion is negligible and the system is essentially 
one-dimensional [14). The dynamics of the BEC is de- 
scribed by the dimcnsionlcss nonlinear Gross-Pitaveskii 
(GP) equation, 

d , fid 2 



+ gW + Kcos65 t (T))<ip, (1) 



where g = 8NaR/r 2 is the scaled strength of nonlinear 
interaction, N is the number of atoms, a is the s-wave 
scattering length, K is the kick strength, 6t(T) represents 
J^S (t — nT), T is the kick period, and 9 denotes the az- 
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imuthal angle. The length and the energy are measured 
in units R and , respectively. The wavefunction 
yj(6,t) has the normalization J Q 27r l^j 2 d9 = 1 and satis- 
fies periodic boundary condition ip (0, t) = ip (0 + 2ir, t). 
Experimentally, the ring-shape potential may be realized 
using two 2D circular "optical billiards" with the lateral 
dimension being confined by two plane optical billiards 
[lfij . The (5-kick may be realized by adding potential 
points along the ring with an off-resonant laser 0. The 
interaction strength g may be adjusted using a magnetic 
field-dependent Feshbach resonance . 

The mean energy of each particle is E(t) = 



ig \ip\ 2 ] ip. To determine the evolu- 
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Figure 1: Plots of average energy E(t) versus the number of 
kicks t for three values of g. The kick strength K — 0.8. 



tion of the energy, we numerically integrate Eq.(l) over a 
time span of 100 kicks, using a split-operator method [r?j |. 
with the initial wavefunction ip being the ground state 
-0 (8, 0) = l/\/27r. The kick period is chosen as T = 2n 
to match the condition for anti-resonance. After each 
kick, the energy E(t) is calculated and plotted as shown 
in Fig.l. We see a remarkable difference among noninter- 
action (Fig. 1(a)), weak interaction (Fig. 1(b)), and strong 
interaction (Fig. 1(c)) cases. For noninteraction, the en- 
ergy E(t) oscillates between two values (anti-resonance) 
and the period is 2T; While in the weak interaction 
case, the amplitude of the oscillation decreases gradu- 
ally to zero and then revives, similar to the phenomena 
of beating in classical waves. The values of the oscilla- 
tion and beat frequencies are obtained by Fourier Trans- 
form and the results are presented in Fig. 2. It is shown 
that, besides the appearance of the beat frequency, the 
interaction also modifies the oscillation frequency; For 
stronger interaction (Fig. 1(c)), i.e. g > 1.96, we hnd that 
the energy's evolution demonstrates an irregular pattern, 
clearly indicating the collapse of the quasiperiodic motion 
and the occurrence of instability. 

For the case of weak interaction, the system is pre- 
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Figure 2: Plots of beat and oscillation frequencies versus the 
interaction strength (a) and kick strength (b), where the scat- 
ters are the results from numerical simulation using GP equa- 
tion and lines from analytic expression Eq.(3). 



dominantly in the two energy levels. We can map the 
system onto a spin model with a two-mode approxima- 
tion to the GP equation [JJJ. By considering the con- 
servation of parity we may write the wavefunction as 
ifi = -j== (a + \plb cos #) , where the populations a and 
b at the ground and excited states satisfy the normaliza- 
tion condition |a| 2 + \b\ 2 = 1. The Hamiltonian in the 
spin representation reads 

n = ~ T + h (5 _ T + ¥ ) + ^ KS * 6 ^ ( 2 ) 

where S z corresponds to the population difference |o| 2 — 
|6| and a,rcta,n(S y / S x ) gives the relative phase a = 
arg(a) — arg(6). This Hamiltonian is similar to a kicked 
top model |lP| , but here the evolution between two kicks 
is more complicated. 

With the spin model, we can readily study the dynam- 
ics of the system. For the case of noninteraction, the evo- 
lution between two consecutive kicks is simply an angle tt 
rotation about the z axis. The spin initially directing to 
north pole stays there for time duration T, then the first 
kick rotates the spin by an angle \[2K about the x axis. 
The following free evolution rotates the spin by an angle 
7r about the z axis. Then, the second kick will drive the 
spin back to north pole through another rotation of y/2K 
about the x axis. With this the spin's motion is two kick 
period recurrence and quantum anti-resonance occurs. 

With interaction, the motion between two consecu- 
tive kicks is approximately described by a rotation of 
7T + g(l + SS Z )/A about the z axis. Compared with the 
noninteraction case, the mean field interaction leads to 
an additional phase shift g(l + 3S Z )/A. This phase shift 
results in a deviation of the spin from S x = plane at 
time 2T~ , i.e., moment just before the second kick. As 
a result, the second kick cannot drive the spin back to 
its initial position and quantum anti-resonance is absent. 
However, the phase shift will be accumulated in future 
evolution and the spin may reach S x = plane at a cer- 
tain time mT~ (beat period) when the total accumulated 
phase shift is ir/2. Then the next kick will drive spin back 
north pole by applying an angle \[2K rotation about the 
x axis. 

The above picture is confirmed by our numerical solu- 
tion of the spin Hamiltonian with fourth-order Runge- 
Kutta method |2fl |. In Fig. 3 we see that, the rela- 
tive phase at the moment just before the even kicks 
increases almost linearly and reaches 27r in a beat 
period. The slope of the increment reads, ^rp = 
(a (4T~) - a (2T~)) /2, which can be deduced analyti- 
cally. With this and through a lengthy deduction, we 
obtain an analytic expression for the beat frequency to 
first order in g, 
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Figure 3: (a) Plots of relative phase versus the number of 
kicks t, where g = 0.1, K = 0.8. (b) Schematic plot of the 
phase shift. nT _ ' + ' represents the moment just before (after) 
the nth kick. 



This expression as well as the relation between the oscilla- 
tion frequency and the beat frequency, f osc = | — ^fbeat, 
is in very good agreement with the numerical results as 
shown in Fig. 2. Therefore the beating provides a method 
to measure interaction strength in an experiment. 

Tuning the interaction strength still larger means en- 
hancing further the nonlinearity of the system. From 
our general understanding of nonlinear systems, we ex- 
pect that the solution will be driven towards chaos, in 
the sense of exponential sensitivity to initial condition 
and random evolution in the temporal domain. The lat- 
ter character has been clearly displayed by the irregular 
pattern of the energy evolution in Fig. 1(c). On the other 
hand, the onset of instability (or chaotic motion) of the 
condensate is accompanied with the rapid proliferation 
of thermal particles. Within the formalism of Castin and 
Dum [lOj , the growth of the number of the noncondensed 
atom will be exponential, similar to the exponential di- 
vergence of nearby trajectories in phase space of classical 
system. The growth rate of the noncondensed atoms is 
similar to the Lyapounov exponent, turning from zero to 
nonzero as instability occurs. 

In Castin and Dum's formalism, the mean number of 
noncondensed atoms at zero temperature is described by 
(SN(t)) = Y^kLi( v k(t)\ v k(t)) , where \v k (t)) are governed 
by 



d 
i— 
dt 



H + gQ\i;\ 2 Q gQ^Q* 
-gQ*r 2 Q -H - gQ*\ip\ 2 Q* 



(4) 

where H = ^ + g\ip\ 2 — M is the chemical potential, 
ijj is the ground state of GP equation and the projection 
operators Q are given by Q — 1 — (ip\. 

We numerically integrate Eq.(4) for the Uk, «fc pairs 
over a time span of 100 kicks, using a split operator 
method, parallel to numerical integration of GP equa- 
tion. The initial conditions \uk(0)), \vk(0)), for initial 
ground state wavefunction ip(9) = 1/V27T, are obtained 
by diagonalizing the linear operator in Eq.(4) |2j|. Af- 
ter each kick the mean number of noncondensed atoms 
is calculated and plotted versus time in Fig. 4(a). We 
find that there exists a critical value for the interaction 
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Figure 4: (a) Semilog plot of the mean number of noncon- 
densed atoms versus the number of kicks t. The thicker lines 
are fitting functions. K = 0.8, g = 0.1 (dashed line, fit- 
ting function O.OOOSi 1 ' 3 ), g = 1.5 (dotted line, fitting func- 
tion 0.001 It 2 ), g — 2.0 (dash dotted line, fitting function 
0.32exp(0.1t)). The inset shows the interaction dependence of 
the growth rate. The scatters are from numerical simulation 
and the solid line is the fitting function 0.33(g — 1.96) 1 / 2 . (b) 
Phase diagram of the transition to instability. 



strength, i.e., g c = 1.96, above which, the mean number 
of noncondensed atoms increases exponentially, indicat- 
ing the instability of BEC. Below the critical point, the 
mean number of noncondensed atoms increases polyno- 
mially. As the nonlinear parameter crosses over the crit- 
ical point, the growth rate turns from zero to nonzero, 
following a square-root law (inset in Fig. 4(a)). This scal- 
ing law may be universal for Bogoliubov excitation as 
confirmed by recent experiments |ll| . 

The critical value of the interaction strength depends 
on the kick strength. For very small kick strength, the 
critical interaction is expected to be large, because the 
ground state of the ring-sha pe B EC with repulsive inter- 
action is dynamically stable [22|. For large kick strength, 
to induce chaos, the interaction strength must be large 
enough to compete with the external kick potential. So, 
in the parameter plane of {g, k), the boundary of insta- 
bility shows a "U" type curve (Fig.4(b)). 

Across the critical point, the density profiles of both 
condensed and noncondensed atoms change dramatically. 
In Fig. 5, we plot the temporal evolution of the density 
distributions of condensed atoms as well as noncondensed 
atoms. In the stable regime, the condensate density os- 
cillates regularly with time and shows clear beating pat- 
tern (Fig. 5(a)), whereas the density of the noncondensed 
atoms grows slowly and shows main peaks around 9 = ±7r 
and 0, besides some small oscillations (Fig. 5(b)). In the 
unstable regime, the temporal oscillation of the conden- 
sate density is irregular (Fig. 5(c)), whereas the density 
of noncondensed atoms grows explosively with the main 
concentration peaks at 9 = ±7r/2 where the gradient den- 
sity of the condensed part is maximum (Fig. 5(d)). More- 
over, our numerical explorations show that the cos 2 9 
mode (Fig5.(b)) dominates the density distribution of the 
noncondensed atoms as the interaction strength is less 
than 1.8. Thereafter, the sin 2 9 mode grows while cos 2 9 
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Figure 5: Plots of condensate and noncondensate densities, 
where K = 0.8. (a,b) g = 0.1; (c,d) g = 2.0. 



mode decays, and finally sin 2 9 mode become dominating 
in the density distribution of noncondensed atoms above 
the transition point (Fig5.(d)). Since the density distri- 
bution can be measured in experiment, this effect can be 
used to identify the transition to instability. 

Finally, although the above discussions have been fo- 
cused on a periodic state of anti-resonance, the transition 
to instability due to strong interactions also follows a sim- 
ilar path for a dynamically localized state jl^l ■ The only 
difference is that we start out with a quasiperiodic rather 
than periodic motion in the absence of interaction. This 
means that it will generally be easier to induce instabil- 
ity but still requires a finite strength of interaction. In 
Fig. 6, we show the nonlinear effect on a dynamically lo- 
calized state at K — 5 and T = 1. For weak interactions 
(.9 = 1) the motion is quasiperiodic with slow growth in 
the number of noncondensed atoms. Strong interaction 
(g = 5) destroys the quasiperiodic motion and leads to 
diffusive growth of energy, accompanied with exponential 
growth of noncondensed atoms that clearly indicates the 
instability of the BEC. Notice that the rate of growth in 
energy is much slower than the classical diffusion rate, 
which means that chaos brought back by interaction in 
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Figure 6: Nonlinear effects on dynamically localized states. 
K = 5, T = 1. (a) Plots of average energy E(t) versus the 
number of kicks t, where dash dotted line corresponds to the 
classical diffusion, g — (dash), g = 1 (dot), g — 5 (solid), 
(b) Semilog plot of the mean number of noncondensed atoms 
versus the number of kicks t. g = 1 (dot), g — 5 (solid). 



this quantum system is still much weaker than pure clas- 
sical chaos. 
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